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ABSTRACT 

Gravitational lensing magnifies the observed flux of galaxies behind the lens. We use 
this effect to constrain the total mass in the cluster Abell 1689 by comparing the lensed 
luminosities of background galaxies with the luminosity function of an undistorted 
field. Under the assumption that these galaxies are a random sample of luminosity 
space, this method is not limited by clustering noise. We use photometric redshift in- 
formation to estimate galaxy distance and intrinsic luminosity. Knowing the redshift 
distribution of the background population allows us to lift the mass/background de- 
generacy common to lensing analysis. In this paper we use 9 filters observed over 12 
hours with the Calar Alto 3.5m telescope to determine the redshifts of 1000 galaxies 
in the field of Abell 1689. Using a complete sample of 146 background galaxies we 
measure the cluster mass profile. We find that the total projected mass interior to 
0.25 /i-i Mpc is M2d{< 0.25 Mpc) = (0.48 ± 0.16) x lO^^^ /j-IMq, where our error 
budget includes uncertainties from the photometric redshift determination, the uncer- 
tainty in the offset calibration and finite sampling. This result is in good agreement 
with that found by number count and shear-based methods and provides a new and 
independent method to determine cluster masses. 

Key vi^ords: galaxies: clusters: general - cosmology: theory - gravitational lensing - 
large scale structure of Universe. 



1 INTRODUCTION 

The use of gravitational lensing as a means of cluster mass 
reconstruction provides a theoretically efficient approach 
without the equilibrium and symmetry assumptions which 
typically accompany virial and X-ray temperature methods. 
Mass determination through application of lens shear proves 
to give good resolution in mass maps although measurement 
of absolute quantities is not possible without external cali- 
bration. This so called sheet-mass degeneracy (Falco, Goren- 
stein & Shapiro 1985) is broken however by methods which 
exploit the property of lens magnification. 

First recognised by Broadhurst, Taylor & Peacock 
(1995, BTP hereafter) as a viable tool for the reconstruction 
of cluster mass, lens magnification has the twofold effect of 
amplifying background source galaxy fluxes as well as their 
geometrical size and separation. This immediately permits 
two separate approaches for measuring lensing mass. The 
first involves selecting a sample of sources with a fiat or 
near-flat number count slope. Magniflcation results in a re- 
duction of their local surface number density owing to the 
dominance of their increased separation over the enhanced 
number detectable due to flux amplification. 



Although contaminated by faint cluster members, Fort, 
Mellier & Dantel-Fort (1997) first reported this dilu- 
tion effect using B and I band observations of the clus- 
ter CL0024-I-1654. Later, Taylor et al. (1998, T98 here- 
after) demonstrated how the dilution in surface number 
density of a colour-selected sample of red galaxies lying 
behind the cluster Abell 1689 enables determination of 
its total mass proflle and 2d distribution. A projected 
mass interior to 0.24 Mpc of M2d[< 0.24 /i"^ Mpc) = 
(0.50 ± 0.09) X 10^^ ft-^Mo was predicted, in good agree- 
ment with the shear analysis of Tyson & Fischer (1995) 
who measured M2d{< 0.24 /i"^ Mpc) = (0.43 ± 0.02) x 
10^^ ^~^M0 and Kaiser (1995) with a measurement of 
M2d{< 0.24/1-1 Mpc) ^ (0.43 ±0.04) x 10^^ h'^MQ. 

Since then, several authors have detected source num- 
ber count depletion due to cluster lensing. Athreya et al 
(1999) observe MS1008-1224 and use photometric redshifts 
to identify a background population of galaxies within which 
they measure depletion. Mayen & Soucail (2000) constrain 
the mass profile of MS1008— 1224 by comparing to simula- 
tions of depletion curves. Gray et al (2000) measure the first 
depletion in the near infra-red due to lensing by Abell 2219. 
Finally and most recently, Rognvaldsson et al. (2000) find 
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depletion in the source counts behind CL0024+1654 in the 
R band and for the first time, in the U band. 

The second mass reconstruction approach permitted by 
magnification forms the primary focus of this paper. The 
amplification of fiux by lens magnification introduces a mea- 
surable shift in the luminosity function of background source 
galaxies. With a sufficiently well defined luminosity func- 
tion derived from an unlensed offset field for comparison, 
this shift can be measured to allow an estimate of the lens 
mass (BTP). This method relies upon a set of observed 
source magnitudes which, if assumed to form an effective 
random sampling of luminosity space, is not limited by noise 
from background source clustering unlike the number count 
method (see Section 5.2 for further discussion). 



This paper presents the first application of mass recon- 
struction using lens flux magnification inferred from the lu- 
minosity function of background samples. Unlike the method 
of T98 who defined their background sample based on colour 
cuts, in this work photometric redshifts of all objects in the 
observed field have been estimated. This not only allows an 
unambiguous background source selection but alleviates the 
need to estimate source distances when scaling convergence 
to real lens mass. 

The following section details the theory of mass recon- 
struction from lens magnification of background source mag- 
nitudes. Section 1^ describes the photometric analysis applied 
to observations of A1689 with the redshifts which result. Ob- 
servations of the offset field which provide the absolute mag- 
nitude distribution required for comparison with the A1689 
background source sample are presented in Section ^ From 
this, a parameterised luminosity function is calculated in 
Section |4.2| necessary for application of the maximum likeli- 
hood method. Following a discussion of sample incomplete- 
ness in Section 5.1, a mass measurement of A1689 is given 
in Section |^ where the effects of sample incompleteness are 
quantified. Finally, a signal to noise study is carried out in 
Section |^ to investigate the effects of shot noise, calibration 
uncertainty of the offset field and photometric redshift error. 



2 MASS RECONSTRUCTION 

Measurement of lens magnification of background source 
fluxes requires a statistical approach in much the same way 
as do shear or number count depletion studies. The basis of 
this statistical method relies on the comparison of the dis- 
tribution of lensed source luminosities with t he lu minosity 
function of an un- lensed offset field. As Section [s^ discusses 
further, for a fair comparison, the population of sources de- 
tected behind the lens must be consistent with the popula- 
tion of objects used to form this un-lensed reference lumi- 
nosity function. 

The effect of lens magnification by a factor on a 
source is to translate its observed magnitude from M to 
M -\- 2.5 logj^g p. In terms of the reference luminosity func- 
tion, <j){M,z), the probability of a background galaxy with 
an absolute magnitude Ai and redshift z being magnified by 
a factor is (BTP) 



<?i(M + 2.51ogio^(z),2) 



(1) 



/ (l,{M + 2.5 log^o^^{z),z)dM' 
Magnification depends on the geometry of the observer-lens 



source system hence for a fixed lens and observer, /i is a 
function of source redshift. This redshift dependence comes 
from the familiar dimensionless lens surface mass density or 
convergence, k{z), and shear, 'y(z), which are related to fj,{z) 
via. 



^^{z) = \[l-^{z)]'-^'{z)\ 



(2) 



We wish to apply maximum likelihood theory using the 
probability in equation (|l]) to determine lens magnification 
and hence k. A parametric luminosity function is therefore 
required and so we take ^(M, z) to be a Schechter function 
(Schechter 1976), 

0(M,z) = 0*(z)lO°-*(*^--*^»^+") exp [_io«'4(«.-M)] . (3) 

The Schechter parameters (f)* , M« and a are determined by 
fitting to the magnitude distribution of the offset field (see 
Section T98 use n as a, likelihood parameter by adopting 
the simplification that all sources lie at the same redshift. 
However this is not possible when each source is attributed 
its own redshift. We must therefore express fj, in terms of a 
redshift depenent quantity and a source-independent likeli- 
hood parameter. 

The most direct solution is to separate the convergence. 
Using the parameter Koo introduced by BTP as the conver- 
gence for sources at z = oo, we can write k,{z) as. 



k{z) = 



Koo/(2), 



y'l + Z-y/l + ZL 



(4) 



We therefore choose Hoo as our likelihood parameter with all 
source redshift dependency being absorbed into the function 
/. The lens surface mass density, E, is then related to Koo 
and the lens redshift, zi,, by 



S(z^) = 



SttG' 



+ 



= 2.75 X 10'"k 



vTTIT- 1. 

(1 + ^.)' 



Zl 



1 



hMeMpc. (5) 



Here, we assume an Einstein-de-Sitter universe for reasons 
of simplicity and because BTP show that this result depends 
only weakly on the chosen cosmological model. 

Before choosing a likelihood function, consideration 
must be given to the shear term in equation (^ . Since shear 
scales with source redshift in the same way as the conver- 
gence, we use the so-called k estimators discussed by T98 
which relate k to 7. At the extremes these are 7 = k for 
the isothermal sphere or 7 = for the sheet-like mass. A 
third variation, motivated by cluster simulations and the 
fact that it has an invertible ^(k) relation, is 7 oc k^''^ (van 
Kampen 1998). This gives rise to the parabolic estimator 
which predicts values of k between those given by the sheet 
and isothermal estimators. Using equation (^, the magnifi- 
cation for these three different cases therefore relates to Hao 



\l-2Ko^f{z)\-^ ISO. 

H{z) = { \[Hirx,f{z) - c] [K,aof{z) - 1 / c] p ^ para. 

[1 — Koo/(z)]^^ sheet 



(6) 



and hence three different estimations of Koo exist for a given 
^. The constant c in the parabolic case is chosen to provide 
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Filter: Ac/AA (nm) 


tint (s) 


Use 


826/137 (1-band) 


6000 


Global SED 


774/13 


6800 


Hq at 2 = 0.185 


703/34 


4100 


Background z 


D14/ZO 




Background z 


572/21 


6300 


Background z 


530/35 


3300 


Background z 


480/10 


4200 


4000Aat z = 0.185 


466/8 


4800 


Ca H,K at z = 0.185 


457/96 (B-band) 


6000 


Global SED 



Table 1. The observations of Abell 1689 in all 9 filters (labelled 
as Ac/AA = central wavelength/FWHM). tint gives the total in- 
tegration time in each filter. The I-band data comes from T98. 



the best fit with the cluster simulations. As in T98, we take 
c = 0.7 throughout this paper. 

The likelihood function for Hoa is then formed from 
equation Aw, 



C[k.oc,) OC ]^P[A/i|^i(«:oo),Zi] 



(7) 



where ^ is one of the three forms in equation (g) and the 
product applies to the galaxies behind the cluster region un- 
der scrutiny. Absolute surface mass densities are then cal- 
culated from Koo using equation (P). 

The probability distribution for n^o obtained from equa- 
tion (|l|) for a single galaxy is typically double-peaked as two 
solutions for Koo exist for a given magnification. The choice 
of peak is determined by image parity such that the peak at 
the higher value of Koo is chosen for a galaxy lying inside the 
critical line and vice versa. The chosen peak is then extrap- 
olated to extend over the full k^d range before contributing 
to the likelihood distribution. In this way, a single-peaked 
likelihood distribution is obtained. 

Evidently, calculation of lens surface mass density in 
this way requires redshift and absolute magnitude data for 
background galaxies together with knowledge of the intrinsic 
distribution of magnitudes from an unlensed offset field. The 
next section details the photometric analysis applied to our 
observations of Abell 1689 to arrive at background object 
redshifts and absolute magnitudes. 



3 PHOTOMETRIC ANALYSIS 
3.1 Data acquisition 

Observations of Abell 1689 were performed with the Calar 
Alto 3.5m telescope in Spain using 8 different filters, chosen 
for photometric distinction between foreground, cluster and 
background objects. In addition, the I-band observations of 
T98 were included to bring the combined exposure time to a 
total of exactly 12 hours worth of useable data characterised 
by a seeing better than 2.1". Table |l| details the summed 
integration time for each filter set together with the motiva- 
tion for inclusion of the filter. Note the narrow band filters 
466/8, 480/10 and 774/13 which were selected to pick out 
spectral features of objects lying at the cluster redshift of 
z = 0.185 (Teague, Carter & Gray 1990). 



Image reduction and photometry was performed us- 
ing the MPIAPHOT (Meisenheimer & Roser 1996) software 
written at the MPIA Heidelberg as an extension to the MI- 
DAS reduction software package. Images were de-biased and 
flattened from typically four or five median filtered dusk sky 
flats observed each night for each filter set. Any large scale 
remnant flux gradients were subsequently removed by flat- 
tening with a second-order polynomial fltted to the image 
background. Cosmic ray removal was carried out using the 
pixel rejection algorithm incorporated in MPIAPHOT. All 
post-reduced images were flattened to a la background flux 
variation of less than 0.02 mag. 



3.2 Galaxy catalogue 

Instead of co-adding images in each fllter set before object 
detection, photometric evaluation was carried out on images 
individually. In this way, an estimate of the uncertainty in 
the photon count for each galaxy could be obtained. The 
mean photon count /f*'™' of a galaxy m observed in a filter 
b was calculated as the usual reciprocal-variance weighted 
sum. 



r(b,m) 



(8) 



where the summation acts over all images belonging to a 
particular filter set and the error on is 



-(6,m) 



-1/2 



(9) 



The quantity (t|'''"' is the standard deviation of background 
pixel values surrounding galaxy m in image i. Background 
pixels were segregated by applying an appropriate cut to the 
histogram of counts in pixels within a box of size 13" x 13" 
(40 X 40 pixels) centred on the galaxy. This cut removed 
the high count pixels belonging to the galaxy itself and any 
other neighbouring galaxies within the box. 

Integrated galaxy photon counts were determined using 
MPIAPHOT which sums together counts in all pixels lying 
inside a fixed aperture of radius 6" centred on each galaxy. 
A 'mark table' accompanying every image in every filter- 
set provided co-ordinates of galaxy centres. Tables were fit 
within an accuracy of < 1" to individual images using copies 
of a master table derived from the deepest co-added image; 
that observed in the I-band. In this way consistent indexing 
of each galaxy was achieved throughout all catalogues. The 
master table was generated using the object detection soft- 
ware 'SExtractor' (Berlin & Arnouts 1996). Only galaxies 
were contained in the master table, star-like objects being 
removed after identification by their high brightness and low 
FWHM. With a detection threshold of 3a above the aver- 
age background fiux, galaxies in the I-band image were cat- 
alogued after coincidence matching with objects detected at 
the 3a level in the associated V-band data presented by T98. 
Despite cataloguing ~ 3000 galaxies, the resulting number 
was limited to a total of ~ 1000 due to the relatively shallow 
data observed with the 466/8 narrow-band filter. 
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3.3 Photometry 



The integration of photon counts in an aperture of fixed 
size requires constant seeing across all images to allow cor- 
rect determination of colours. To ensure constant seeing, all 
images were degraded by Gaussian convolution to the worst 
seeing of 2.1" measured in the 466/8 filter before galaxy 
counts were evaluated. The effects of changing weather con- 
ditions were compensated for by normalising images within 
each filter set to an arbitrarily chosen image in that set. 
Normalisation was conducted by scaling the galaxy counts 
in each image so that the average counts of the same stars 
in all images was equal. This ensured correct calculation of 
the weighted counts and the error from equation (^) and 
(^). These quantities were later scaled to their calibrated 
photometric values. 

Calibration of the photometric fluxes from the weighted 
counts was provided using the spectrum of the dominant 
elliptical galaxy in the centre of A1689 taken from Pickles & 
van der Kruit (1991). Denoting this spectrum as the function 
Fs{X), the calibration scale factors kb for each filter b were 
calculated using 



r(6,s) 



kb / dA 



E{\)Ti,(\)Fs{\)\ 
he 



(10) 



where the function Tb(A) describes the filter transmission 
efficiency, E[\) is the combined filter-independent efficiency 
of the detector and telescope optics and 7'^'''''' is the mea- 
sured integrated photon count rate of the central galaxy. 
The values of kb obtained in this way were only relatively 
correct owing to the lack of an absolute calibration of the 
published spectrum. Absolute calibration scale factors were 
calculated from observations of the standard star G60-54 
(Oke 1990) in all filters in exactly the same manner. Verifi- 
cation of this absolute calibration was provided by the con- 
sistency of ratios of fc;, (absolute) /fci, (relative) to a zero point 
of Am = 2.11 ± 0.01 magnitudes averaged over all filters. 

Consideration of equation (^) shows that only the 
quantity 



dA£;(A)rf,(A)F„(A)A = 



(11) 



can be known for any galaxy m with a calibrated photon 
count rate. The required photometric flux 



F 



{b,m) 



dAri,(A)F„(A) 



(12) 



can not therefore be directly determined without making 
an approximation such as the simpliflcation of filter trans- 
mission curves to top hat functions. Although this is ac- 
ceptably accurate in narrow band filters, it is not for broad 
band filters. This problem was avoided by the more sophisti- 
cated technique of fitting model spectra to measured galaxy 
colours as the next section discusses. 



3.4 Photometric redshift evaluation 

Direct calculation of photometric fluxes using equation ( p^ 
was made possible by fitting model spectra to the set of cal- 
ibrated photon count rates measured for each galaxy across 
all filters. Expressed more quantitatively, equation ( pT| ) was 
applied for each filter to a library of template spectra to 




Figure 1. Redshift distribution of the 958 objects photometri- 
cally evaluated in the field of A1689. The grey histogram plots 
all 958 redshifts whereas the black histogram plots only the 470 
redshifts with a la error in redshift of less than 0.05. The peak 
at z ~ 0.18 is the contribution from the cluster galaxies. 



arrive at a set of scaled fllter counts for each spectrum. 
Galaxies were then allocated library spectra by flnding the 
set of library colours which best fit the measured galaxy 
colours. Note that this differs from conventional template 
fitting where spectra are redshifted and scaled to fit observed 
colours in a much more time costly manner. 

The spectral library was formed from the template 
galaxy spectra of Kinney et al. (1996). A regular grid of 
galaxy templates was generated, varying in redshift along 
one axis from z — to z = 1.6 in steps of Az — 0.002 and 
ranging over 100 spectral types from ellipticals, through spi- 
rals to starbursts along the other. 

The set of photometric errors given by equation (^ for 
an individual galaxy across all filters gives rise to an error el- 
lipsoid in colour space. Using the size and location of these 
error ellipsoids, probabilities of each library entry causing 
the observed sets of colours for each galaxy were then cal- 
culated as 



p{q\2 



v/(2^)"l^ 



cxp 



1 

2 ^ 



[q,-Q,{z,s)f 



(13) 



where n is the number of colours, ctj comes from propagation 



of the error given by equation (|9|) and V = diag(cri, cr^). 
Each galaxy's position vector in colour space, q = (gi, g„) 
is compared with the colour vector Q of the library spec- 
trum with a given redshift z and type s. Finding the maxi- 
mum probability corresponding to the closest set of match- 
ing colours therefore immediately establishes redshift and 
galaxy type. An assessment of the uncertainty in this red- 
shift is subsequently obtained directly from the distribu- 
tion of the probabilities associated with neighbouring library 
spectra. 

Figure |^ shows the distribution of the 958 successfully 
classified galaxy redshifts estimated from the full filter set. 
We measure an average redshift error of (a^) = 0.08. The 
maximum redshift limit of z ~ 0.8 comes from the condition 
that the 4000A limit must lie in or blue-ward of the second 
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reddest filter in tiie set. Tlie pealc at 2 ~ 0.18 in Figure ^ is 
clearly the contribution from the cluster galaxies. 

The feature at z ~ 0.4 is most likely real and not 
an artifact of the photometric method. Such artifacts oc- 
cur due to 'redshift focusing' when particular redshifts are 
measured more accurately than others. Where the uncer- 
tainty is larger, galaxies can be randomly scattered out of 
redshift bins, producing under-densities and corresponding 
over-densities where the redshift measurement is more ac- 
curate. This effect depends on the details of the filter set, 
being more common when fewer filters are used, but can be 
modelled by Monte Carlo methods. 

The top half of Figure ^ shows the results of one re- 
alisation of such a Monte Carlo test for redshift focusing. 
The plot indicates how accurately the method reproduces 
redshifts of spectra scaled to I = 20 with photometric noise 
levels taken from the A1689 filter set. Each point repre- 
sents a single library spectrum. Reproduced spectral red- 
shifts, Zphot, were determined by calculating colours through 
application of equation to the library spectra with red- 
shifts 2iib. These colours were then randomly scattered by an 
amount determined from the filter-specific photometric error 
measured in the A1689 data before application of the red- 
shift estimation method outlined above. The bottom half of 
Figure ^ shows the same plot generated using spectra scaled 
to I = 21 with the same photometric error taken from the 
A1689 data. 

The accuracy of reproduced redshifts at I = 20 is clearly 
better than those at I = 21 where photometric noise is more 
dominant. The lack of any sign of redshift focusing in the 
vicinity of 2; ~ 0.4 leads us to conclude that the feature 
seen at this redshift in the A1689 data is probably real. 
The I = 21 plot which corresponds approximately to our 
sample magnitude cut of B = 23.6 (see Section 5.1) shows 
that galaxies at redshifs below z = 0.05 have on average 
higher estimated redshifts. This only marginally affects the 
overall redshift distribution and yet partly explains the lack 
of galaxies at z < 0.05 in the A 1689 redshifts of Figure |l|. It 
is worth emphasising here that the significance of the peak 
at « ~ 0.18 attributed to the cluster galaxies is far in excess 
of any effects of redshift focusing. 

Figure ^ shows a comparison of the photometrically de- 
termined redshifts Zphot around the peak of the redshift dis- 
tribution of Abell 1689, with spectroscopically determined 



(Teague, Carter & Gray, 1990) redshifts 



A very slight 



bias between Zphot and Zspcc can be seen. This bias is quanti- 
fied by fitting the line Zphot = Zs-pcc+c by least-squares to the 
data points which gives c = 0.0036. Referring to equation 
( psj ) shows that if this small bias is applied to all redshifts 
in our sample, a negligible difference of Akoo = 0.001 would 
result. Our filter set was selected primarily to distinguish the 
cluster members, hence at higher redshift we must rely on 
our Monte Carlo estimates of the redshift uncertainty (see 
Section 6.1). 

Abell 1689 lies in a region of sky where there is a very 
low level of galactic dust. Our redshifts are therefore not 
affected by this source of contamination. However, dust in 
the cluster itself is another concern. We have modelled the 
effects of reddening by cluster dust and find that although 
magnitudes are slightly affected, the redshifts remain the 




• •' iiiiiiPi 

■iHlli 

ill I i: !: i 



Ill ii 
li 



II 

I*. : 



liN'^:!::;.. 
ill I- 



0.2 



0.4 



0.6 



0.8 



Figure 2. A single Monte Carlo realisation showing the accu- 
racy of the photometric redshift evaluation method. Input library 
spectra with redshifts are scaled to I = 20 (top) and I = 21 
(bottom) and subsequently used to calculate sets of colours us- 
ing the A1689 filterset. These colours are randomly scattered by 
the filter-specific photometric errors measured in the A1689 data 
before calculating the reproduced redshifts Zphot- 



4 OFFSET FIELD AND LENS CALIBRATION 

The unlensed, intrinsic magnitude distribution required by 
the likelihood analysis was taken from an offset field ob- 
served as part of the Calar Alto Deep Imaging Survey 
(CADIS) conducted by the Max-Planck Institut fiir As- 
tronomie, Heidelberg (Meisenheimer et al. 1998). Data for 
this survey were observed to a complete depth of B ~ 24.5 
mag in 16 filters from the B-band to the K-band with the 
2.2m telescope at Calar Alto. We use the CADIS 16-hour 
field for our mass calibration. 

Using exactly the same methods outlined in the previ- 
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0.18 



0.2 



0.22 



Figure 3. Comparison of the photometric redshifts estimated 
in the cluster Abell 1689 with spectroscopically determined red- 
shifts. The distribution shows slight non-Gaussianity in the error 
distribution. The mean redshift of the cluster determined spec- 
troscopically is z = 0.185 (Teague, Carter & Gray 1990), while 
the mean photometric redshift is z = 0.189 ± 0.005. 



ous section for the A1689 data, photometric redshifts and 
rest-frame absolute magnitudes were determined for all ob- 
jects in the field. In addition to galaxy templates in the 
spectral library however, quasar spectra from Francis et al. 
(1991) and stellar spectra from Gunn & Stryker (1983) were 
also included, the primary motivation for this difference 
being the CADIS quasar study. As a by-product, a more 
sophisticated object classification method was achieved by 
finding the overall object class yielding the highest signifi- 
cant probability given by equation ([l^). Details of this and 
the CADIS quasar study are given in Wolf et al. (1999). To 
ensure a fair comparison between the ofi'set field and the 
cluster field, the CADIS B-band (Ac/AA = 461/113 nm) 
galaxy magnitudes were used with the A1689 B-band mag- 
nitudes in the likelihood analysis discussed in Section ^. 

Investigation of evolution (see for example, Lilly et al. 
1995, EUis et al. 1996) of the CADIS luminosity function 
is left for future work. A preliminary study indicated no 
significant evolution which would impact on the lens mass 
determination. We therefore applied the same redshift se- 
lection as used for the Abell 1689 background sources, and 
assumed a no-evolution model. 

We present two estimates of the calibration B-band lu- 
minosity function: a nonparametric 1/Knax method (Section 
4.1) and a maxim um likelihood parametric fit to a Schechter 
function (Section 



4.2 



The former method allows us to see 
the distribution of luminosities without imposing a precon- 
ceived function, and gives a visual impression of the uncer- 
tainties in the parametric fit. In addition the Knax approach 
allows us to make bas ic t ests for sample completeness (we 
discuss this in Section 5.1). The latter maximum likelihood 



fit provides a convenient function for performing the sec- 
ond likelihood analysis to determine lens magnification. We 
begin by describing the nonparametric method. 



4.1 The nonparametric CADIS B-band 
luminosity function 

An estimate of the luminosity function of galaxies in the 
CADIS B band was provided initially using the canonical 
1/Knax method of Schmidt (1968). The quantity Knax is 
computed for each galaxy as the comoving volume within 
which the galaxy could lie and still remain in the redshift 
and magnitude limits of the survey. For an Einstein-de-Sitter 
universe, this volume is. 



V -(^ 

\ -no 



5uj 



min(2u ,z-r, 
max(z; ,Zm, 



dz 



(1 + 2)3/2 



(14) 



where Suj is the solid angle of the observed field of view and 

D{z) is 

i3(2) = 2(l-(l + 2)-i/^). (15) 

The upper limit of the integral in equation ( |l^ is set by the 
minimum of the upper limit of the redshift interval chosen, 
Zu, and the redshift at which the galaxy would have to lie to 
have an apparent magnitude of the faint limit of the survey, 
-2m„-,ax- Similarly, the maximum of the lower limit of the 
chosen redshift interval, zi, and the redshift at which the 
galaxy would have to lie to have an apparent magnitude of 
the bright limit of the survey, 2m„i„ , forms the lower limit of 
the integral. This lower integral limit plays a non-crucial role 
when integrating over large volumes originating close to the 
observer where the volume element makes only a relatively 
small contribution to Knax. 

The redshifts 2m„ax s-nd 2m„,i„ are calculated for each 
object by finding the roots of 



M ~ miini 4- 51ogi(, 



{l + z)D{z) 



ho 



-K{z) = -42.39 (16) 



where M is the absolute magnitude of the object, mum is the 
appropriate maximum or minimum survey limit and K{z) 
is the K-correction. Although the K-correction for each ob- 
ject at its actual redshift was known from its apparent and 
absolute rest-frame magnitude, the redshift dependence of 
this K-correction was not. In principle, this redshift depen- 
dence could have been calculated directly for each object 
using its best-fit spectrum returned from the photometric 
analysis. However, the much simplified approach of approx- 
imating the K-correction as a linear function of redshift was 
employed. This was primarily motivated by its improved ef- 
ficiency and the relatively weak influence the K-correction 
was found to have on the final luminosity function. Lilly 
et al. (1995) find that the K-correction in their B-band for 
elliptical, spiral and irregular galaxies is proportional to red- 
shift in the redshift range < z < 0.8. The following form 
for K{z) was thus adopted 



K{z) 



Kjzo) 
zo 



(17) 



where K{zo) is the K-correction of the object at its actual 
redshift zo- 

Once Vmax has been calculated for all objects, the maxi- 
mum likelihood estimated luminosity function </!> at the rest- 
frame absolute magnitude M in bins of width dM is then. 



(t){M)dM = 



1 



(18) 
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Figure 4. The CADIS B band object luminosity function calcu- 
lated with the l/Vmax formalism. Errors account only for errors 
in redshift. Points lie at bin centres, the widths of each chosen to 
hold the same number of objects. There are 371 galaxies in total 
selected by the redshift limits 0.3 < z < 0.8 and the apparent B 
magnitude me < 24. 5. T he solid line is the Schechter function 
determined in Section 4.2 



where the sum acts over all objects with magnitudes between 
M - AM/2 and M + AM/2. 

Figure ^ shows the luminosity function of B band mag- 
nitudes from the CADIS offset field which has a solid view- 
ing angle of Su) = 100 arcmin^. To match the selection of 
objects lying behind A1689, only objects within the redshift 
range 0.3 < z < 0.8 were chosen. A further restriction on 
the apparent B magnitude of me < 2 4.5 w as applied for 
completeness of the sample (see Section 5.1), yielding a to- 
tal of 371 galaxies. The data points in Figure W are centred 
on bins chosen to maintain an equal number of objects in 
each. 

The Icr errors shown here were calculated from Monte 
Carlo simulations. Object redshifts were randomly scattered 
in accordance with their associated errors provided in the 
CADIS dataset. For each realisation, the Knax of each ob- 
ject was re-calculated using the re-sampled redshift. The 
resulting standard deviation of the distribution of values of 
for each bin given by equation (|l^) was then taken as 
the error. In this particular instance, no consideration was 
given to the magnitude errors or the propagation of the red- 
shift error into object magnitudes. Section W/A discusses this 
further. 



4.2 Parameterisation of the CADIS B-band 
luminosity function 

The maximum likelihood method of Sandage, Tammann & 
Yahil (1979, STY hereafter) was employed to determine the 
Schechter function best describing the CADIS B-band mag- 
nitudes. This parameterisation is essential for the determi- 
nation of lens mass using the likelihood approach. 

In much the same way as the probability in equation 
(^) was formed, the STY method forms the probability pi 
that a galaxy i has an absolute magnitude Mi, 



Pi = p{Mi\zi) 



min(Mn: 
max(M^ 



e,x(zi),M2) 
.n(^,).Mi) 



(f){M)AM 



(19) 



where Mmax(zi) and Afmin(zi) are the absolute magnitude 
limits corresponding to the apparent magnitude limits of 
the survey at a redshift of Zi. Conversion of these apparent 
inagnitude limits includes the K-correction using equation 
(td) with z set to Zi. A further restriction is placed upon 
the integration range by imposing another set of magnitude 
limits Ml < M < M2 which for the CADIS data were set 
at the maximum and minimum absolute magnitudes found 
in the sample. 

The likelihood distribution in this case is a two di- 
mensional function of the Schechter parameters M« and a 
formed from the product of all probabilities pi. The best fit 
M* and a are therefore found by maximizing the likelihood 
function. 



\n£.{PL,a) = 
Eti{ln<^(M.; 



In/ 



min(Mn: 
max(Mn, 



.x(z,),M2) 
,(zi),Mi) ' 



■AM} +Cp 



(20) 



with the constant Cp arising from the proportionality in 
equation (p^). An estimate of the errors on M* and a are 
calculated by finding the contour in a, M* space which en- 
compasses values of a and M* lying within a desired confi- 
dence level about the maximum likelihood /Imax. 

To account for uncertainties in the Schechter parame- 
ters due to the redshift and magnitude errors derived by the 
photometric analysis, Monte Carlo simulations were once 
again performed. Redshifts and absolute magnitudes of the 
entire sample were randomly scattered before re-application 
of the maximum likelihood process each time. The error in 
absolute magnitude was calculated using simple error prop- 
agation through equation ( p^ yielding. 



lU 5_ 1-0.5(1+^0 

Zi In 10 1 + — (1 + 2, 



-1/2 
W2 



2 2 
0"zi + C^m 



(21) 



for each object with redshift Zi and apparent magnitude m^. 
Here, the K-correction given by equation (^) has been used 
such that the quantity Ki = K{zi) is calculated from equa- 
tion (0) using rrii, Mi and zi as they appeared in the CADIS 
dataset. The typical ratio of apparent magnitude to redshift 
error was found to be crmi/o'zi ~ 1% due to the relatively 
imprecise nature of photometric redshift determination. 

The final errors on M, and a were taken from an ef- 
fectively convolved likelihood distribution obtained by com- 
bining the scattered distributions produced from the Monte 
Carlo simulations. Figure ^ shows the final la and 2a like- 
lihood contours calculated allowing for redshift and magni- 
tude errors. These predict the resulting parameters 



M, = -19.43 



+0.47 

0.64 



+ 5l0gi(,/l. 



1 ^r-l-0.25 

a = -l.45_o.23 



(22) 



where projected la errors are quoted. For completeness, the 
normalisation was also calculated to be 



= 0.0164/1^ Mpc"^ 



(23) 



in this case where evolution of the luminosity function has 
been neglected. 
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Apparent Magnitude 



Figure 5. Likelihood contours for tlie CADIS B-band Schochter 
parameters taking photometric redshift and magnitude error into 
consideration. Icr and 2a confidence levels are plotted correspond- 
ing to Aln/3 = 1.15 and 3.09 respectively. 



5 SAMPLE CONSISTENCY 

Ensuring that the A1689 sample of sources is consistent with 
the CADIS offset field sample is necessary to prevent biases 
from entering our results. The first level of compatibility 
we have already enforced by applying a redshift selection of 
0.3 < z < 0.8 to both samples. The second, discussed in 
Section 



5.1 



below, is sample completeness. A slightly less 
obvious consideration must also be given to galaxy morpho- 
logical type as Section 5.2 explains. 




Apparent Magnitude 



5.1 Completeness 

Determination of the faint magnitude limit beyond which 
both the A1689 and CADIS data set become incomplete 
is important for the calculation of an accurate lens mass. 
Both samples must be complete for fair comparison. Incor- 
rect evaluation of the CADIS limiting magnitude results in 
larger values of Knax and hence a biased luminosity function 
not representative of the intrinsic A1689 distribution. Simi- 
larly, completeness of the A1689 sample also affects the lens 
mass in a manner quantified in Section 



6.5 



Estimation of the completeness of both data sets was 
provided using the l^/Vmax statistic (Schmidt 1968). In this 
ratio, Vmax is calculated using equation (^4|) whereas V is 
the comoving volume described by the observer's field of 
view from the same lower redshift limit in the integral of 
equation ( p^ ) to the redshift of the object. If a sample of 
objects is unclustered, exhibits no evolution and is complete, 
the position of each object in its associated volume Knax will 
be completely random. In this case, the distribution of the 
V/Vmax statistic over the range to 1 will be uniform with 
(V/V^s.,) = 0.5. 

If the sample is affected by evolution such that more 
intrinsically bright objects lie at the outer edges of the Vmax 
volume, then V^/Vmax is biased towards values larger than 
0.5. The reverse is true if a larger number of brighter ob- 
jects lie nearby. If the sample is incomplete at the limiting 
apparent magnitude chosen, estimations of Vmax will be on 
average too large and will cause V/Vmax to be biased towards 



Figure 6. Variation of y/Vmax with limiting apparent B magni- 
tude for the A1689 (top) and the CADIS (bottom) sample. The 
grey region corresponds to the la errors described by equation 
1241) which are an underestimate due to the unconsidered effects 
of galaxy clustering. 



values less than 0.5. The requirement that (V/Knax) =0.5 
for completeness is also subject to fluctuations due to finite 
numbers of objects. In the absence of clustering, the un- 
certainty due to shot noise on (V/Knax) calculated from A'' 
galaxies can be simply shown to be 



12A' 



(24) 



In order to arrive at an apparent magnitude limit for 
the CADIS and A1689 fields, values of (V/Vmax) were calcu- 
lated for different applied limiting magnitudes and plotted as 
shown in Figure The grey region in this plot corresponds 
to the la errors described by equation (|2^) which lessen at 
the fainter limiting magnitudes due to the inclusion of more 
objects. Clustering adds extra noise and so these errors are 
an underestimate of the true uncertainty (one can show that 
the uncertainty in V/Vmax increases by the square root of 
the average number of objects per cluster). 

Without knowledge of the effects of clustering. Fig- 
ure ^ shows that a limiting magnitude of ttib < 24.5 for 
CADIS and tub < 23.6 for A1689 corresponds to a value 
of (V/Vmax) — 0.5 and thus completeness. These magni- 
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Figure 7. The Abell 1689 (top) and CADIS 16'" (bottom) field 
redshift-apparent B magnitude parameter space. The dashed box 
highlights the selection criteria for the 146 A1689 background 
galaxies. 



A1689 to contain a higher proportion of ellipticals than a 
sample of field galaxies away from a cluster environment. 
Since the luminosity function of E/SO galaxies is thought 
to be different from that of spirals (see, for example, Chiba 
& Yoshii 1999 and references therein), comparison of our 
A1689 sample with the CADIS offset field sample may be 
expected to introduce a bias. 

A related issue stems from the fact that the determina- 
tion of an object's photometric redshift requires its detection 
in every filter belonging to the filter set. Both the CADIS 
and A1689 filter sets contain narrow band filters which cause 
the main restriction on which objects enter into the photo- 
metric analysis. Our test in the B-band therefore 
does not give a true limiting B magnitude but one which 
applies only to objects detectable across all filters. As long 
as both samples under comparison are complete according 
to this test, the sole consequence of this is that certain mor- 
phological types will be under-represented. With the CADIS 
filter set of 16 filters differing from the 9 filters used for the 
observation of A1689, this again might be expected to cause 
inconsistent galaxy types between fields, biasing results. 

Fortunately, our photometric analysis yields morpho- 
logical type in addition to redshift. We are therefore able to 
directly measure the fraction of galaxies of a given type in 
both samples and thus test for biases. We find that in the 
CADIS sample, the ratio of E/SO:Spiral galaxies is 12%±7% 
and for the A1689 sample, this is 25% ± 13% with Pois- 
son errors quoted. This is in reasonable agreement with the 
canonical Postman and Geller (1984) E/SO:S fraction for 
field galaxies of about 30%. Approximately 60% — 70% of 
both samples are classified as starburst galaxies. Given the 
uncertainty in these fractions, we would argue that as far as 
we can tell, they are consistent with each other. Without a 
bigger sample of galaxies and possibly spectroscopically con- 
firmed morphologies we are unable to do better although as 
the measurements stand, we would not expect any serious 
inconsistencies in morphology which would bias our results. 



tudes are in agreement with the apparent magnitude limits 
at which the number counts begin to fall beneath that mea- 
sured by deeper surveys (such as the B band observations 
of Lilly, Cowie & Gardner (1991) which extend to a depth 
of TTiB — 26) and also correspond to a IOct object detection 
threshold deduced from the photometry. 

Figure Q shows the distribution of galaxies in the 
redshift-apparent magnitude plane in the Abell 1689 field. 
Abell 1689 itself can clearly been seen at z = 0.18. Superim- 
posed is the region of parameter space we use for the mass 
determination with 0.3 < z < 0.8 and 18 < niB < 23.6 
yielding 146 galaxies. 

5.2 Morphological Type 

Perhaps the most difficult inconsistency to quantify is that of 
variation in galaxy morphological type between samples. It 
has been known for some time that elliptical galaxies cluster 
more strongly than spirals (Davis & Geller 1976) and that 
the elliptical fraction in clusters is an increasing function 
of local density (Dressier 1980). One might therefore expect 
a sample of galaxies lying behind a large cluster such as 



6 MASS DETERMINATION 

Taking the CADIS luminosity function as a good estimate of 
the intrinsic distribution of A1689 source magnitudes in the 
range 0.3 < z < 0.8, we can calculate and E assuming 
a relation between k and 7. As noted in Section ^, in the 
absence of a detection of evolution of the luminosity func- 
tion with redshift, we assume a no-evolution model for the 
background sources in this field. In general, the occurrence 
of evolution is anticipated however we expect its inclusion 
in the background model to have only a minor effect on the 
derived mass. 

For the purpose of comparison with other studies we 
shall quote values of k as well as Koo and E. As k is de- 
pendent on the source redshift, this is not a useful quantity 
to quote when the redshift distribution is known. The con- 
vergence we quote is the redshift averaged quantity defined 
by 

(25) 

i = l 

where Nt, is the number of source galaxies. For the field 
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of Abell 1689 we find that (/b) 
source redshift of z^a = 0.45. 



0.57, giving an effective 



6.1 Sources of uncertainty 

Three sources of error on were taken into consideration: 

1) The maximum hkelihood error obtained from the width 
of the Ukehhood distribution at In/Imax — 0.5Ax^, with Ax'^ 
the desired confidence leveL AU object magnitudes and red- 
shifts were taken as presented directly in the A1689 data 
while assuming the Schechter parameters from equation 

2) The uncertainty of the Schechter parameters i\f, and a 
from the likelihood analysis of the CADIS offset field. 

3) The redshift and magnitude uncertainties of individual 
objects in the A1689 data, derived from the photometric 
analysis. 

In Section ^ we will show that the contribution of each source 
of uncertainty to the overall error depends on the number 
of galaxies included in the analysis. Taking all 146 galaxies 
across the entire field of view, the errors from each con- 
tributor listed above, expressed as a percentage of the total 
standard deviation were found to be; 50% from the maxi- 
mum likelihood (essentially the shot noise), 25% from the 
uncertainty in M» and a and 25% from the redshift and 
magnitude error. 

The latter two sources of error in the above list were si- 
multaneously included using the Monte Carlo method. 1000 
simulations were carried out, randomly drawing values of M* 
and a from the convolved likelihood distribution shown in 
Figure ^. For each realisation, redshifts and absolute magni- 
tudes of objects in the A1689 field were scattered in exactly 
the same fashion as before with the CADIS dataset using 
their associated photometric errors. The standard deviation 
of the scattered values of produced in this way was then 
added in quadrature to the uncertainty of the maximum 
likelihood error obtained from item one of the list above to 
give the overall error. 

The magnitude calibration error of gam ~ 0.01 dis- 
cussed in Section B.3 was ignored. Inspection of the form of 
the Schechter function in equation shows that a system- 
atic magnitude offset is exactly equivalent to an error in M, . 
Clearly, the la error quoted for M« in equation ( |22| ) com- 
pletely overwhelms this magnitude calibration uncertainty 
which was therefore deemed insignificant. 

Finally, the dependence of our measurement of itoo on 
the feature seen in the A1689 redshift distribution at 2 ~ 0.4 
was tested. We removed all galaxies contributing to this peak 
and re-calculated the results of the following two sections. 
Apart from a larger uncertainty due to the decreased num- 
ber of objects, we found very little difference from the results 
obtained from the full dataset, indicating that our measure- 
ment is not dominated by the concentration of galaxies at 
z ~ 0.4. 



6.2 The differential radial k profile 

Background source objects from the A1689 data set were 
binned in concentric annuli about the cluster centre for the 



calculation of a radial mass profile. The relatively small num- 
ber of objects contained in the sample however was unfortu- 
nately insufficient to allow computation of a profile similar 
in resolution to that of T98. 

Apart from the effects of shot noise, this limitation re- 
sults from the simple fact that bins which are too narrow 
do not typically contain a large enough number of intrin- 
sically bright objects. This has the effect that the knee of 
the Schechter function assumed in the likelihood analysis is 
poorly constrained. As equation (^) shows, a large uncer- 
tainty in M, directly results in a large error on the magni- 
fication and hence on Koo ■ Experimentation with a range of 
bin widths quickly showed that in order to achieve a tolera- 
ble precision for Hoo, bins had to be at least ~ 1.1 arcmin in 
width. With the observed field of view, this gave a limiting 
number of three bins, illustrated spatially in the lower half 
of Figure ^. In Section 5.c we find the average profile within 
an aperture, which provides a more robust measurement of 

K. 

The top half of Figure ^ shows the k data points. These 
were converted from the maximum likelihood derived Koc for 
each bin using equation (|2^). The profile of T98 is shown 
superimposed for comparison. Upper points correspond to 
the sheet estimator while the lower points are due to the 
isothermal estimator. The Icr error bars plotted wer e ca lcu- 
lated taking all three contributions listed in Section 6.1 into 
account . 

Despite relatively large errors, the data points show an 
amplitude in good agreement with the profile derived from 
the number count study of T98. These errors seem large 
in comparison to those of the number count profile but the 
number count errors do not take the systematic uncertainties 
in background count normalisation, number count slope or 
background source redshift into consideration. 

It is noticeable that the data points suggest a profile 
that is perhaps a little fiatter than that derived by T98. It 
appears that more mass is detected at larger radii although 
this is not particularly significant. 



6.3 Aperture k profile 

In addition to the radial profile, the variation of average 
surface mass density contained within a given radius can be 
calculated. By applying the likelihood analysis to the ob- 
jects contained within an aperture of varying size, a larger 
signal to noise can be attained at larger radii where more 
objects are encompassed. With a small aperture, the same 
low galaxy count problem is encountered as Figure ^ shows 
by the large uncertainty in this vicinity. In this plot, the 
parabolic estimator of equation is used to obtain Koo 
since as T98 show, this is a good average of the isothermal 
and sheet estimators and agrees well their self-consistent axi- 
symmetric solution. Application of the axi-symmetric solu- 
tion is not viable in this case since we are limited to only 3 
bins. 

Using equation (|2^), Koo is again scaled to k. The grey 
shaded region in Figure |^ depicts the la errors, with the 



sources of uncertainty from Section 6.1 taken into account. 
The thin solid and dashed black lines show the variation 
of aperture k and its error calculated by averaging the 
parabolic estimator profile presented in T98. The error does 
not account for uncertainties arising from background count 
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Figure 8. Top Comparison of radial k profiles for A1689. Data 
points show isothermal (lower) and sheet (upper) estimated k 
obtained from this work, la error bars are plotted. The solid 
lines indicate the same isothermal/sheet estimator-bound profile 
obtained by T98 using integrated number counts with Icr errors 
shown by dashed lines. Bottom: Spatial location of the annular 
bins on the A1689 field of view. Open dots are objects selected 
by 2 > 0.2 and solid dots by z < 0.2. 



normalisation, number count slope or background source 
redshift. The results of T98 were shown to be in good agree- 
ment with the shear analysis of Kaiser (1996) and hence we 
find a consistent picture of the mass amplitude and slope 
from all three independent methods. 



As expected from the results of Section 3.2, generally 
more mass than that predicted from the number counts is 
seen, especially at large radii. The following section quanti- 
fies this for a comparison with the projected mass result of 
T98. 



6.4 Projected mass 

From t he v alues of n^o used to generate the «: profile in 
Section |6.2| and the result of equation the cumulative 
projected masses in Table were calculated. Errors were 



Figure 9. Variation of average surface mass density contained 
within a given radius R (thick dark line). Icr errors are shown by 
the grey shaded region. The thin black solid and dashed lines show 
the average surface mass density and Icr errors (due to likelihood 
analysis only) of the parabolic estimated profile of T98. 



Radius (arcsoc) 


M2d(< R) 


65 


(0.16 ±0.09) X IO^/i-^Mq 


130 


(0.48 ±0.16) X lQ^"'h-^MQ 


195 


(1.03 ±0.27) X lO^^fe-iMe 



Tab le 2 . Cumulative projected mass given by the profile of Sec- 
tion \>'\ 

derived from propagation of the errors on the binned values 

of Koo. 

These projected masses are in excellent agreement 
with those of T98. At the redshift of the cluster 1' = 
Q.IVJ Mpc and hence the second cumulative mass listed 
in Table |2| gives 

M2d{< 0.25 Mpc) = (0.48 ± 0.16) x lO^'^ H'^Mq (26) 

which is perfectly consistent with the result from the number 
count study. The error here is comparable to the 30% error 
quoted for the result of T98 when allowing for all sources of 
uncertainty. The projected mass contained within 195 arc- 
sec is a little higher than that predicted by T98 although re- 
mains arguably consistent given the errors involved in each. 



6.5 Effects of sample incompleteness 

One final uncertainty not taken into consideration so far is 
that of sample incompleteness. Changing the limiting ap- 
parent B magnitude in the determination of the CADIS lu- 
minosity function directly affects the fitted values of A/*, 
a and hence the maximum likelihood Kqo. Similarly, differ- 
ing numbers of objects included in the A1689 sample from 
variations in its limiting magnitude also has an influence on 

Table |^ quantifies this effect for the CADIS objects. It 
can be seen that increasing the faint limit mmax (ie. includ- 
ing fainter objects) has little effect on Koo until the limit 
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'^max 


A/, 


a 




25.5 


-19.06 


-0.80 


0.61 


25.0 


-19.25 


-1.10 


0.65 


24.5 


-19.43 


-1.45 


0.70 


24.0 


-19.98 


-1.87 


0.74 


23.5 


-19.24 


-1.53 


0.75 



(iso) Koo(para) Kqo (sheet) 



+0.03 
-0.04 
+0.04 
-0.04 
+0.06 
-0.04 
+0.03 
-0.04 
+0.04 
-0.04 



3+0.04 
-0.06 
7+0.06 
-0.06 
-+0.08 
-0.08 
91+013 

90+" 0'' 



0.69: 
0.77j 
0.85^ 



-,+0.07 
-0.08 
^+0.08 
-0.09 
K+0.10 
-0.10 
lOg+0.16 

1.10 



0.76: 
0.84^ 
0.96"^ 



-0.17 
+0.10 
-0.09 



Table 3. Variation of limiting apparent B magnitude rrimax of the 
CADIS field and its effect on the Schechter parameters and the 
resulting value of Koo • Values of M» assume h = 1 . The apparent 
magnitude limit of ft = 23.6 was assumed for the A1689 data in 
calculating the maximum likelihood k,oo- Errors here are taken 
only from the width of the likelihood curves. 



• •• •• .r:-.-\'^;' ' 

' *• •* • * • ^ , • , • • ' 

f.v.--*! »;'■'?«•.'•. 





Koo(iso) 


too (para) 


too (sheet) 


24.5 
24.0 
23.5 


r, ^7 + 0.03 
'-0.03 

'^•"J-0.04 

U.OM_o.07 


1 03+"" 
i^-^'^-O.ii 

94+" i0 

79+0.10 
'^-0.11 


1 30+015 
l-^'J-0.13 
1 -19+0.13 

'J-y^-0.12 



Table 4. Variation of the maximum likelihood determined Kq 
with limiting apparent B magnitude m^, ^ of the A1689 data. 



The Schechter parameters of Section 4.2 were assumed in the 
likelihood analysis. 



rriniax — 24.5 is reached. Beyond this limit, Koo starts to fall. 
Two inferences can therefore be made . Fi rstly, this suggests 
that the magnitude limit in Section 5.1 from the V/Vmax 
test, being consistent with the limit here, was correctly cho- 
sen. Secondly, tiao is relatively insensitive to the choice of 
fWmax if the sample is complete (and not smaller than the 
limit at which shot noise starts to take effect). 

The effect of varying the magnitude limit of the A1689 
sample is quantified in Table ^ A clear trend is also seen 
here; as its limiting faint magnitude mmax is reduced, Koo 
falls. Assuming linearity, a rough estimate of the uncertainty 
of Kaa given the uncertainty of the sample magnitude limit 
is given by: 



{~ O.lAmmax isothermal 
~ 0.2Ammax parabolic (27) 
~ 0.4Ammax sheet 

Referring to Figure |^, a suitable uncertainty in mmax of 
the A1689 sample of say ±0.2 magnitudes might be argued. 
If this were the case, the projected masses of the previous 
section calculated with the parabolic estimator would have 
a further error of ~ 5% which is negligibly small. 



7 SIGNAL-TO-NOISE PREDICTIONS 

Including all possible contributions of uncertainty in the cal- 
culation of mass, the previous section showed that even with 
relatively few galaxies, a significant cluster mass profile can 
be detected. One can make predictions of the sensitivity of 
the method with difi'ering input parameters potentially ob- 
tained by future measurements. This exercise also serves as 
an optimisation study, enabling identification of quantities 
requiring more careful measurement and those which play 
an insignificant part. 



Figure 10. Correlation of photometric redshift error with appar- 
ent B-band magnitude for the A1689 data. No significant corre- 
lation between cr^ and z exists. 



The most convenient means of carrying out this inves- 
tigation is by application of the reconstruction method to 
simulated galaxy catalogues. Catalogues were therefore con- 
structed by randomly sampling absolute magnitudes from 
the Schechter function fitted to the CADIS offset field in 
Redshifts were assigned to each magnitude by 



Section 4.2 



randomly sampling the distribution parameterised by T98 
(their equation 22) from the Canada France Redshift Sur- 
vey (Lilly et al. 1995). A range of catalogues were pro- 
duced, varying by the number of objects they contained and 
their distribution of galaxy redshift errors modeled from the 
A1689 data. 

Figure ^ shows how the distribution of photometric 
redshift error, a^, correlates with the A1689 B-band appar- 
ent magnitude. No significant correlation between cr^ and 
redshift was found. Catalogue objects were thus randomly 
assigned redshift errors in accordance with their apparent 
magnitude, given by the correlated distribution in Figure 
[l^ . Different catalogues were generated from different scal- 
ings of this distribution along the CTz axis. 

Each catalogue was then lensed with a sheet mass char- 
acterised by Koo = 1 before applying the reconstruction. 
1000 Monte Carlo realisations were performed for each re- 
construction, scattering object redshifts according to their 
assigned errors in the same manner as in the reconstruction 
of A1689. Furthermore, to model the uncertainty associated 
with the offset field, assumed values of the Schechter pa- 
rameters M, and a were once again subject to Monte Carlo 
realisations. All catalogues were reconstructed assuming sets 
of Schechter parameters drawn from a range of scaled ver- 
sions of the distribution shown in Figure ^ 

The resulting scatter measured in the reconstructed 
value of Koo for each catalogue and assumed a-M« scaling 
was combined with the average maximum likelihood error 
across all realisations of that catalogue to give an overall 
error. This total error was found to be well described by, 

2 ^ 1 + (2^^^) + (0.12aM. f + (0.37a„)' - 0.18cr<,M. (28) 
n 

where n is the number of galaxies, cr^ is the sample average 
redshift error and gm, and aa are the projected errors on M» 
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and a respectively as quoted in equation (^2|). The quantity 
faAf. is the covariance of a and M* defined by 



/ C{M,,a)[M,- {M-,))[a-{a)) dM.da 



(29) 



/£(M.,a) dA/.da 

where the likelihood distribution C is given by equation ( po| ) . 
We find that Gcm, = 0.039 for the CADIS offset field. Equa- 
tion (||) is valid for n > 20 and 0.0 >g^> 0.3. 

Equation ( p8| ) shows that when the number of objects 
is low, shot noise dominates. With n ~ 200 however, un- 
certainties from the calibration of the offset field start to 
become dominant. The factor of 2 in the photometric red- 
shift error term stems from the fact that redshift errors also 
translate directly into absolute magnitude errors through 
equation (^). Another discrepancy arises when comparing 
this redshift error with the redshift error contribution of 
25% claimed for the A1689 data in Section 3.1. This is ac- 



counted for by the fact that K-corrections were present in 
the A1689 data whereas in the simulated catalogues there 
were not. Equation jil] ) quantifies the increase in magnitude 
error with the inclusion of K-corrections. This translates to 
an approximate increase of 20% in the overall error with an 
average K-correction of —1.0 for the Af689 data. 

Emphasis should be placed on the criteria for which 
equation (|2^) is valid. The predicted overall error rises dra- 
matically when fewer than ~ 20 objects are included in the 
analysis. Simulations with 15 objects resulted in maximum 
likelihood errors rising to beyond twice that predicted by 
simple shot noise. Th is stems mainly from the effect men- 
tioned in Section 6.2, namely the failure of the likelihood 
method when the knee of the Schechter function is poorly 
constrained. 

The most immediate improvement to a multi-colour 
study such as this would therefore be to increase galaxy 
numbers. As previously noted, only when bins contain ~ 200 
objects do offset field uncertainties become important. Ob- 
serving in broader filters is one way to c ombat the limit pre- 
sented by galaxy numbers. Section 3.2 noted how the 3000 



galaxies detected in the 1 band image were instantly reduced 
to 1000 by the shallow depth limit placed by the narrow 
466/8 filter, even though both were observed to similar inte- 
gration times. Using broader filters will also inevitably give 
rise to less accurate photometric redshifts. However as the 
analysis of this section has shown, one can afford to sacrifice 
redshift accuracy quite considerably before its contribution 
becomes comparable to that of shot noise. 

Deeper observations provide another obvious means of 
increasing the number of galaxies. The error predictions 
above indicate that the expected increase in galaxy num- 
bers using an 8 metre class telescope with the same expo- 
sure times as those in this work should reduce shot noise by 
a factor of ~ 3. Since deeper observations would also reduce 
redshift and offset field calibration uncertainties to negligi- 
ble levels, the only source of error would be shot noise. In 
this case, the signal to noise for from equation (^8|) be- 
comes simply Koo \/n and hence our mass estimate for A1689 
could be quoted with a 9(t certainty. 



8 SUMMARY 

Photometric redshifts and magnitudes have been deter- 
mined for objects in the field of Abell 1689 from multi- 
waveband observations. This has allowed calculation of the 
luminosity function of source galaxies lying behind the clus- 
ter. Comparison of this with the luminosity function ob- 
tained from a similar selection of objects in an unlensed off- 
set field has resulted in the detection of a bias in the A1689 
background object magnitudes attributed to lens magnifica- 
tion by the cluster. 

To ensure that systematic biases do not affect our re- 
sults, we have given careful consideration to the consistency 
between both the A1689 dataset and the CADIS offset field 
dataset. We find the distribution of galaxy types within the 
redshift range 0.3 < z < 0.8 applied to both samples to be 
very similar. This demonstrates that our lower redshift cut- 
off is sufficient to prevent objects in the A1689 dataset from 
being significantly influenced by the cluster environment. 

After allowing for sources of uncertainty due to redshift 
and magnitude error, offset field calibration error and like- 
lihood error (including shot noise), a significant radial mass 
profile for A1689 has been calculated. We predict a projected 
mass interior to 0.25 Mpc of 

M2d{< 0.25 h-^ Mpc) = (0.48 ± 0.16) x 10^^ /i'^Mq (30) 

in excellent agreement with the number count analysis of 
T98 and the shear results of Tyson & Fischer (1995) and 
Kaiser (1995). 

We can compare the efficiency of the method presented 
in this paper in terms of telescope time and signal-to-noise 
with the number count method used by T98. The 5.5(j re- 
sult quoted by T98 does not include uncertainty due to 
their background count normalisation, number count slope 
or background source redshift distribution. Adding these er- 
rors to their result gives an estimated signal-to-noise of 3(j, 
the same as this work. Regarding telescope time, we define 
a 'total information content' for each study as the prod- 
uct of telescope collecting area and total integration time. 
T98 observed 6000s in each of the V and I bands with the 
NTT 3.6m. Comparing with the CA 3.5m telescope and 12 
hours integration time used in this study shows that we 
have amassed a total information content of approximately 
3 times that required for the T98 result. 

Despite the extra time penalty induced by our method, 
we note that deeper observations, especially in the narrow 
band filters used, would increase the signal-to-noise of our 
result significantly since we are dominated by shot noise. 
Our signal-to-noise analyses in Section |^ showed that a 9(t 
detection of mass would be possible using an 8 metre class 
telescope, equivalent to an increase in integration time by a 
factor of 5. This is in contrast to T98 whose main source of 
uncertainty comes from the unknown source redshift distri- 
bution. Shot-noise makes a negligible contribution to their 
error to the extent that increasing their total integration 
time by a factor of three to match the total information 
content of this work would still result in a signal-to-noise of 
3(7 (Dye 1999). 

This paper has been primarily devoted to establishing 
the viability of lens mass reconstruction using the luminosity 
function method. We have shown that the two main advan- 
tages over the number count method employed by T98 are 
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that use of photometric redshifts have enabled breaking of 
the mass/background degeneracy and that the technique is 
independent of their clustering if it is assumed that they 
form an effective random sampling of luminosity space. 
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